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The retrieval of phases from intensity measurements is a key process in many fields in science, 
from optical microscopy to x-ray crystallography. Here we study phase retrieval of a one-dimenslonal 
multi-phase object that is illuminated by quantum states of light. We generalize the iterative 
Gerchberg-Saxton algorithm to photon correlation measurements on the output plane, rather than 
the standard intensity measurements. We report a numerical comparison of classical and quantum 
phase retrieval of a small one-dimensional object of discrete phases from its far-held diffraction. 
While the classical algorithm was ambiguous and often converged to wrong solutions, quantum light 
produced a unique reconstruction with smaller errors and faster convergence. We attribute these 
improvements to a larger Hilbert space that constrains the algorithm. 


I. INTRODUCTION 


Quantum states of light have been widely explored 
in recent years for their ability to offer considerable en¬ 
hancement in measurement sensitivity over classical ones 
[l|. Quantum states were mostly considered for enhanc¬ 
ing the sensitivity of the measurement of a single op¬ 
tical phase using an optical interferometer. Recently, 
the problem of simultaneous estimation of several ot- 
tical phases using quantum light was investigated Q . 
Here we investigate an iterative phase-retrieval technique 
for the estimation of one-dimensional phase objects with 
quantum states of light, and show that the technique is 
more robust than its classical version. 

The problem of phase retrieval is one of great scien- 
tihc interest, which arises when the intensity recorded in 
the far-held is used to determine the phase structure of 
an object, information which is otherwise undetected. It 
stems from the fact that detectors record the intensity of 
waves, while often the important information is encoded 
in their phases. Phase retrieval has been intensively in¬ 
vestigated 0 and found applications across many helds 
of science, from astronomy (wave-front sensing) Q to 
nanotechnology (x-ray crystallography and electron mi¬ 
croscopy) 0|- In the most common scenario, the far- 
held diffraction intensity pattern is measured; This mea¬ 
sured far-held intensity, together with known constrains 
on the illumination (e.g. its intensity prohle) are used in 
an iterative algorithm to derive the phase structure of the 
object Q. It is known, however, that phase retrieval has 
certain limitation; for example, phase-retrieval of one¬ 
dimensional objects is problematic, and often leads to 
multiple ambiguous solutions 

Here we study the use of quantum states of light to 
measure a one-dimensional multi-phase object, as shown 
schematically in Fig. [TJ using a phase-retrieval approach. 
We present a protocol that utilizes entangled states to 
reconstruct multiple phases simultaneously using an it¬ 
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erative error-reduction algorithm. It is shown by our 
numerical results that a quantum approach has a few 
advantages over the approaches that use classical light. 
First, phase retrieval using quantum light can be unam¬ 
biguous, as it reaches the single, correct solution. The 
algorithm also converges much faster than in the clas¬ 
sical case. Furthermore, the use of quantum states can 
enhance the sensitivity of the retrieved phases over cases 
where classical light is used, for the same number of pho¬ 
tons probing the system. These quantum enhanced ca¬ 
pabilities are already revealed for two-photon entangled 
states, and are particularly important when probing deli¬ 
cate samples which are sensitive to illumination intensity, 
such as biological samples quantum gases , and 
atomic ensembles [ll|. With recent advances in gener¬ 
ating quantum states in x-ray [H, their application for 
phase retrieval holds great promise. 

An important milestone in the field of phase retrieval 
was the iterative algorithm of Gerchberg and Saxton Q , 
that reduces the error in the phase object with every iter¬ 
ation. It is doing so by iterating between the input plane 
Ein{x), and the output plane Eout{u), related by the 
known transformation of the system, and applying the in¬ 
formation obtained from the intensity of the input plane 
and the output intensity measurement at each iteration 



FIG. 1: A schematic description of phase retrieval using quan¬ 
tum states. A quantum state |i/)), entangled over m modes, 
is input to a multi-mode interferometric system. The state 
passes through m phases denoted hy 6 = { 61 , 62 ,... , 0 m,}, 
and followed by a transformation U. Photon correlation mea¬ 
surement is carried out on the output state. 
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of the algorithm. In the common case of far-held diffrac¬ 
tion, the transformation between the input and output 
planes is the Fourier transform, where x and u are the co¬ 
ordinates along the input and output planes, respectively. 
Similarly, phase retrieval can be applied to the frequency¬ 
time domains, where the temporal held Emit) is related 
to its spectral one Eoutiyj) by a Fourier transformation, 
and the problem relates to determining temporal am¬ 
plitudes and phases from spectral power measurements 
[Flj . However, phase reconstruction is not unique in gen¬ 
eral. As the amplitudes of the input and output planes 
\Ein{x)\, and \Eout{u)\ are restricted by their measured 
intensities, there might be additional solutions for the 
phase image which are incorrect. Some examples include 
shifted images Einix — a:o), mirror images E*^{—x), and 
global phases e^^Ein{x), which are exact solutions of the 
phase retrieval problem, yet there may exist many other 
non-trivial exact solutions as well 14, [i3 , which to the 
best of our knowledge were never analyzed. 

Quantum enhanced phase retrieval uses quantum 
states to probe the object, as we outline in the next sec¬ 
tion. In the far-held, instead of intensity measurements, 
we employ measurement of photon correlations on the 
output quantum states. Using these measurements and 
the knowledge of the input state, we describe the algo¬ 
rithm used for the retrieval of the phases of the object. 
We then describe a specihc example, where we also dis¬ 
cuss in some detail the sensitivity of using quantum light 
for phase retrieval. 


II. METHODOLOGY 

Let us consider the problem of estimating a one¬ 
dimensional object of multiple phases, probed by a quan¬ 
tum state of light, as shown in Fig. [T] A phase ob¬ 
ject is characterized by a set of m unknown phases, 
9 = {9i,... ,9m}- 


A. Quantum light 


operator for mode x. The state in Eq. [ijthen becomes 

D 

( 2 ) 

fc=i 

where the set of phases accrued by the state, (j) = {(j)k} 
is related to the object phases 9 hy 4>k = 9 ■ . 

Next, this quantum state undergoes a transformation, 
most commonly Fourier transformation via diffraction, 
which transforms the state in Eq. [5] to the final state at 
the output 


D 

IV'f) = (3) 

i=l 

In Eq. 13] we assumed that the transformation described 
by U is unitary; it can therefore be represented by a 
unitary m x m matrix U. Using such operation trans¬ 
forms the photon creation operators, in any mode x, by 
dl —>■ One example of such an operation 

is the discrete Fourier transform (DFT), which we use 
as the transformation to the far-held plane. The DFT is 
represented by an m x m matrix \U\x,y = exp(i27r(x — 
l)(2/-l)/m)/Vm [Hi. The set of amplitudes denoted by 
P in Eq. |3| can be calculated as a function of the input 
state amplitudes a of Eq. lU 


D 


Pt = = E 


afce^’^'°Per(Vfc^t) 




( 4 ) 


where 14 is an N x N sub-matrix of the matrix U con- 

(k) 

structed by repeating the xth row of C/ rix times, and 
then repeating the yth column of the result matrix riy ^ 
times for all x and y, and Per(I4_t) is the permanent of 
the matrix Vk^t- 

Finally, all D probabilities in the output state Pp = 
{Pp}\ = {|/3tP} are measured by employing Wphoton 
coincidence detection in m modes, as described in Fig. 

[II 


An initial pure state of N photons in m modes has the 
form 

H) = ^ak\r}^\nf \ ... ( 1 ) 

k^l k^l 

where is a vector of length m with photon number 
components n'x'^ in each mode x and for each conhg- 
uration k, such that = N. The set of am¬ 

plitudes d = {offc} (where k = 1,...,D) is normalized 
HfeLi = 1, and the total number of configurations 

After passing through the phase object the state ac¬ 
crues m phases, as described by the unitary transforma¬ 
tion Ug- = exp(^ 02,713,), where fix is the number 


B. Phase retrieval algorithm 

We will now describe the procedure for reconstruct¬ 
ing the set of unknown phases 9, which generalizes the 
Gerchberg-Saxton (GS) error reduction iterative algo¬ 
rithm Q to quantum light, as shown in FiglDJa). The goal 
is to retrieve the phases 9 from the known amplitudes of 
the input state a, and the measured probabilities of the 
output state Pp. 

We begin by guessing random initial values for the 
set of m phases, The zth iteration of the algo¬ 
rithm begins by constructing an input state = 

as in Eq. |21 where = 

^(*-1) . nW, Then, the state is transformed as 
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FIG. 2: A schematic description of phase retrieval algorithm 
for (a) quantum and (b) classical light. Each iteration in the 
algorithm uses the input state amplitudes, a or Etn, trans¬ 
forms these amplitudes {U), and applies the measured pho¬ 
ton correlations Pp or intensities lout, which is followed by 
the inverse transformation (17^), for the quantum or classical 
algorithms respectively. The object phases 9 evolve over the 
iterations of the algorithm, while for the quantum algorithm 
these phases are found from 4 >- 

in Eq. [3] to find a set of output state amplitudes 
The arguments of these complex amplitudes are 
combined with the measured output probabilities Pp 
to yield a new estimate of the output state \ip^p) = 
®xp(^arg(/3(*^))|f^(^)). This state is then 
transformed back to retrieve the corresponding input 
state = J2k=i exp(*(/)^*^)|n('=)), from 

which a new estimate for the set of phases 0*-*^ is found 
by inverting the relation 

The GS algorithm is known to always converge to a so¬ 
lution by means of error reduction [8|, llTj , not necessarily 
to the correct solution. In order to quantify the error 
with which the algorithm converges we use two different 
measures: the error in the Fourier output state in the ith 
iteration, SPpii), 

and the phase error in the ith iteration, S9{i)^ 

m 

= mod (|6»W - 6»,,|,27r)2, (6) 

x=2 

where mod (a, 6) is the reminder of the division a/b 
rounded to the nearest value. 


C. Required conditions for uniqueness 

In order to achieve a unique phase retrieval the input 
quantum state amplitudes d of Eq. [T] are chosen such 
that they satisfy two conditions. 


1. Avoiding trivial ambiguities 

In order to eliminate the trivial ambiguities, the in¬ 
put states should be chosen such that they have no sym¬ 
metries of translation and reflection with respect to the 
phases of the modes. The simplest solution is to arrange 
the average photon number in these modes to break those 
symmetries. 

2. Phase transformation 

The object phases 9 are evaluated by first estimating 
the phases (j) of the quantum states that are used as the 
interrogating field. To uniquely determine the m phases, 
clearly one has to start with at least m basis states am¬ 
plitudes d that are non-zero. When m initial amplitudes 
are used, the matrix that expresses the relations between 
the m object phase 9 and the subset of m phases that are 
used from (j) should have non-zero determinant, so that 
the phases of each mode can be uniquely extracted from 
the reconstructed phases of the basis vectors. In addi¬ 
tion, since the phases of the state 4> are reconstructed up 
to an integer number of 27r, we need to make sure that 
this shift remains an integer number of 27r for the object 
phases 9 as well. 

III. EXAMPLE: RETRIEVAL OF SIX PHASES 

We describe here a phase retrieval problem with m = 6 
as an instructive example. We assume that the transfor¬ 
mation U is the discrete Fourier transform (DFT) which 
is the most relevant one for many practical realizations. 
We begin with the following quantum two-photon state 
(A = 2): 

l^fe) = ^( 12 , 0 , 0 , 0 , 0 , 0 ) + 11 , 1 , 0 , 0 , 0 , 0 ) + 

| 1 , 0 , 1 , 0 , 0 , 0 ) + | 1 , 0 , 0 , 1 , 0 , 0 ) + 

|1,0,0,0,1,0) + |0,1,0,0,0,1)). (7) 

The input state of Eq. [7] contains only six of the 
Dm=6 = 21 configurations of two photons in six modes, 
all six with equal amplitudes of 1/ VE. This input state 
was constructed with care in order to fulfil the two re¬ 
quirements for uniqueness 1section fllCl) : First, the state 
is chosen to eliminate trivial ambiguities, i.e. translation 
or reflection. Indeed, the input state of Eq. [7] has non¬ 
equal average photon numbers in the various ports; the 
intensity ratio between the modes is6:2:l:l:l:l, 
which breaks both symmetries. Second, the input state 
basis is chosen such that it enables the extraction of mea¬ 
sured object phases 9, from the phases of the quantum 
state (j), which are actually estimated by the algorithm. 

To compare the phase retrieval performance with 
quantum light to that with classical light, we per¬ 
formed two sets of simulations, with the same set of 
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object phases that was chosen randomly to be 9obj = 
{0,3.22,4.10,4.57,1.35,4.11}. 



FIG. 3: Different solutions found by the phase retrieval algo¬ 
rithm with classical light. Solutions (a)-(g) are wrong recon¬ 
structions, while (h) is the correct one. 
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FIG. 4: Comparing the performance of the phase retrieval 
algorithm with qnantum and classical light, (a) Histograms 
of the retrieved phase error SO and (b) the Fourier errors SPf 
as a function of the iteration number i, using classical and 
quantum light for retrieval of Oobj for 1000 runs of the algo¬ 
rithm. In the classical case, only 16% correct reconstructions 
were achieved, while the erroneous solutions are the majority 
of the instances. The Fourier error of the classical algorithm 
is shown only for cases that converged to the correct 
solution. The quantum case, that used the entangled two- 
photon state given in Eq. [T] always converged to the correct 
phases. 

First, we applied the GS algorithm with classical light, 
as shown in Fig[2}b). We assumed a classical coherent in¬ 
put field with input amplitudes Ein = {VG, V^, 1,1,1,!}, 
such that it reproduces the same intensity ratio as the 
quantum state in Eq. [3 and, similarly, does not have 
reflection and translation symmetries. This input field 
was transformed with the phases 0obj and then by the 
DFT to obtain a set of six complex output amplitudes 

Eout- The output intensities lout are then used as the 
input to the GS algorithm. We ran the algorithm a large 
number of times, each run starting with a different ran¬ 
dom set of initial phases. The algorithm almost always 
converged, i.e. found a solution which reproduces the in¬ 
tensities in the Fourier plane with very low error in the 
Fourier plane, SPp ^ 10“^, but most of the times it 
did not find the correct set of phases Oobj- All the solu¬ 
tions that were found using classical light are presented 
in Fig. [31 where 7 out of 8 of these solutions are actually 
wrong. A histogram showing the phase error distribution 


for 1000 runs of the algorithm is shown in Fig. HJa). In 
this representative example, the algorithm converged to 
the correct phases only in about 16% of the runs. 

In contrast, the quantum algorithm always found the 
correct solution. The quantum simulation was performed 
in an analogous way: the input state of Eq. [7] was first 
transformed by the phase vector Oobj according to Eq. 
13 and then by the DFT to calculate the set of output 
amplitudes /3, as in Eq. 01 Note that in contrast to 
the m = 6 amplitudes that characterized the classical 
case, here there are 0^=6 = 21 amplitudes, that describe 
all the combinations of two photons in six modes. The 
values of .P /3 = are used as input to the GS algorithm, 
that was run many times with random initial phases. As 
shown in Fig. |4}a), it always converged to the correct 
phase vector Oobj ■ 

Even when the algorithm using the classical light con¬ 
verged to the correct phases (i.e. in about 16% of the 
times), it did so less efficiently than the quantum one. 
Fig. mb) shows the progressive reduction of Fourier 
plane error SPp with increasing iterations for both the 
quantum and classical states of light, averaged over many 
runs. For the quantum state, this error is given by Eq. 
m and similarly, for the classical state it is given by 
(^^(c0)2 ^ ^ ^here 

Ex is the far-held amplitude of the xth mode in the ith 
iteration El- 


IV. SENSITIVITY 

In the simulations described above, the far-held am¬ 
plitudes, either classical or quantum, were calculated 
from theory. In practice, these values would be mea¬ 
sured. Measurements with classical or quantum states of 
light are often limited in sensitivity due to shot-noise. In 
fact, metrology with nonclassical states of light is more 
often than not motivated by its superior sensitivity in 
phase measurements. While this is not the focus of this 
work, for completeness, we wish to compare the preci¬ 
sion of the algorithms with classical and quantum light. 
For this purpose, we performed Monte-Garlo simulations 
with quantum and classical light, for the same example 
that we discussed in the previous section, given a to¬ 
tal number of photons passing through the sample Nt- 
Again, in the classical case, we considered only the runs 
that yielded the correct phases which were less than 16% 
of the total runs. In practice, of course, there is no way to 
identify the correct solution, but here we are interested 
to check the ultimate precision of the algorithm. 

The phase error 50 for the quantum and classical light 
are shown as a function of Nt in Fig. |S] Also shown 
in the figure is the minimal error that could be achieved 
with classical light, = (jn — 1)/\/Nt 0. The 

phase error 56, achieved with the particular quantum 
state of Eq. |7| and for the particular phase object Oobj, 
can be fitted by 8.0/-/Nt- From Fig. |5]it is evident that 
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FIG. 5: Phase error for the algorithm using quantum and 
classical light input 59, and the ultimate minimal error limit 
achievable with classical light 59^°'^^^, as a function of total 
number of photons probing the system, Nt- Here, too. The 
classical line is drawn only for the subset of runs, 16% of 
all runs, that converged to the correct solution. 

this phase error of the quantum state is better than that 
achieved with classical light, 11.8/\/Nt, again pointing 
to the relative advantage of the quantum algorithm. 

Our particular state does not show sensitivity below 
classical limits (the ultimate classical limit in our exam¬ 
ple is b/y/Nr). Still it performs 48% better than classical 
light with the same intensity distribution, again, even if 
compared only with the small fraction of classical solu¬ 
tions that converged to the correct solution. We note 
that a method for calculating the ultimate limit of sensi¬ 
tivity for any input pure state is given in Ref. Q. Fur¬ 
thermore, enhancement factors of sensitivity greater than 
was shown here are probably possible, either by consid¬ 
ering higher photon numbers N > 2 then in our example, 
or by optimizing the input state amplitudes, as well as 
the transformation U. 


V. DISCUSSION 

The solution found with quantum light is unambigu¬ 
ous, converges faster, and is more precise than the one 
found using classical light. The reason for this lies in 
the fact that the quantum states are characterized by 
D = amplitudes, which is significantly larger 

than the m amplitudes of the classical case for any N > 1, 
21 vs. 6 in our two-photon in six modes example above. 
This significant increase in the number of constraints im¬ 
posed on the Fourier plane quantum state amplitudes 
for the same number of unknown phases is most proba¬ 
bly what leads to the elimination of extra solutions, the 
faster convergence, and the very small error. 


We have considered in the example presented above the 
phase retrieval of 6 modes, however, we have checked that 
our approach performs perfectly also for higher number 
of modes. We tested our quantum algorithm with two 
photon states {N = 2) for m = 10, 20, and 30 unknown 
phases. For that, we have generalized the quantum state 
of Eq. [7] by introducing additional terms of the form 
|1,0,..., 0,1, 0,..., 0) having in total m terms of equal 
amplitudes, still meeting the symmetry requirements for 
the input state. The algorithm performed just as well, 
retrieving the phases accurately for all values of m tested, 
while the classical algorithm with coherent light with in¬ 
tensities that matched that of the generalized quantum 
state had much reduced rate of successful reconstruction, 
less than 1% for m > 10. 

Finally, we discuss the considerations for practical re¬ 
alization of the suggested method. First, arbitrary entan¬ 
gled states of light are generally hard to generate, how¬ 
ever, probabilistic generation of multi-mode correlated 
two-photon states is technically possible, although, to the 
best of our knowledge, has never been demonstrated in 
multi-mode systems, perhaps due to lack of interest. The 
DFT can be implemented by Fourier multiport devices, 
which have been ejmerimentally demonstrated for quan¬ 
tum states of light [3, [3 • Additionally, measurements of 
fV-photon correlations over m modes requires a set of D 
measurements which for large number of modes/photons 
can be challenging. A similar problem has been recently 
encountered in the realizations of the boson-sampling 
problem [Tol - I^ . For a quantum state of two photons 
however, the problem involves only m{m + l)/2 mea¬ 
surements, which is quite practical, using, for example, 
large arrays of single photon detectors |23j . or cameras 
with single-photon sensitivities . 

It is important to note here that quantum measure¬ 
ment of correlations on classical coherent input light will 
not be useful: they will not yield any additional informa¬ 
tion or any other advantage, as classical states are uncor¬ 
related and separable, unlike quantum states which can 
exhibit inherent photon correlations between the modes. 
We also note that these findings raise many interesting 
theoretical questions, for example on the optimal choice 
of the input quantum state and on its ultimate sensitivity, 
and how the technique will perform with two-dimensional 
objects, as well as objects of phase and absorption. 


VI. SUMMARY 

We studied the use of quantum states of light for phase 
retrieval. We showed that quantum states of two photons 
exhibit a few advantages in retrieving the phase of a one¬ 
dimensional object from its far-held diffracted intensity, 
as compared with classical states of light. The quan¬ 
tum approach achieves a unique reconstruction, which 
converges faster, and is more robust when subjected to 
shot-noise, when compared with classical approaches. 
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